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Context. The ubiquitous presence of the Fe line complex in the X-ray spectra of galaxy clusters offers the possibility of measuring 
■^L ■ their redshift without resorting to spectroscopic follow-up observations. In practice, the blind search of the Fe line in X-ray spectra is 

a difficult task and is affected not only by limited S/N (particularly at high redshift), but also by several systematic errors, associated 
with varying Fe abundance values, ICM temperature gradients, and instrumental characteristics. 

Aims. We assess the accuracy with which the redshift of galaxy clusters can be recovered from an X-ray spectral analysis of Chandra 
archival data. We present a strategy to compile large surveys of clusters whose identification and redshift measurement are both based 
on X-ray data alone. 

' Methods. We apply a blind search for K-shell and L-shell Fe line complexes in X-ray cluster spectra using Chandra archival obser- 

vations of galaxy clusters. The Fe line can be detected in the ICM spectra by simply analyzing the C-statistics variation ACatat as a 
function of the redshift parameter, when all the other model parameters are frozen to the best-fit values. We repeat the measurement 
under different conditions, and compare the X-ray derived redshift zx with the one obtained by means of optical spectroscopy Zo. 
^ ' We explore how a number of priors on metallicity and luminosity can be effectively used to reduce catastrophic errors. The ACstat 

I provides the most effective means for discarding wrong redshift measurements and estimating the actual error in zx ■ 

, Results. We identify a simple and efficient procedure for optimally measuring the redshifts from the X-ray spectral analysis of clus- 

■ ters of galaxies. When this procedure is applied to mock catalogs extracted from high sensitivity, wide-area cluster surveys, such as 

those proposed with Wide Field X-ray Telescope (WFXT) mission, it is possible to obtain complete samples of X-ray clusters with 
reliable redshift measurements, thus avoiding time-consuming optical spectroscopic observations. Our analysis shows that, in the case 
, of WFXT, a blind Fe line search is 95% successful for spectra with more than 1000 net counts, whenever ACstat > 9, corresponding 

^ ■ formally to a 3a confidence level. The average en'or in the redshift zx decreases rapidly for higher values of ACatat- Finally, we 

0^ ' discuss how to estimate the completeness of a large cluster samples with measured zx- This methodology will make it possible to 

00 ' trace cosmic growth by studying the evolution of the cluster mass function directly using X-ray data. 

O ! 

^ Key words, galaxies: clusters: general - cosmology: observations - X-ray: galaxies: clusters - intracluster medium 

O • 

1 . Introduction regions of groups and clusters, has proven to be very efficiency 

in investigating the chemical properties of the ICM. Indeed, the 

The study of clusters of galaxies allows one to investigate the ubiquitous presence of the K-shell Fe line complex at 6.7-6.9 

_ _ large-scale structure of the Universe, constrain the cosmological ^.^y^ ^^^-^^^ ^^^^ ^^^^^^^^ by Mitchell et al. (1976), has been 

parameters and the spectrum of the primordial density fluctua- ^^^^^^^^ ^^^^^^^ ^^^^^^^^ z ^ 1.4 where X-ray clusters 

^ ; tions, and study the interactions between the member galaxies ^^^^ ^een observed (see Rosati et al. 2004; Stanford et al. 2005; 

, and the ambient intra cluster medium (ICM, see Rosati et al. r^^^^j 2009). In principle, abundance measurements based 

2002; Schuecker 2005; Voit 2005; Borgani 2006; Tozzi 2007; ^^e detection of the Fe line also provide a way to measure the 

Vikhhnm et al. 2009; Mantz et al. 2010; Bohnnger & Wemer ^^^^1^^^ ^he Fe line, thus the cluster redshift. 

2010). The X-ray band is an optimal observational window ^ f . . -f. u u j • v » 1 

. , . ■', , . , ^ r , ■ r^, , , A few redshifts have been measured using X-ray spectral 

for identifying and studying clusters of galaxies. The thermal , . f ■ .. , » • u 

/, ° . . , ° , , T^. ... analyses in cases of no previous optical spectroscopic observa- 

bremsstrahlung emission due to the hot diffuse ICM is roughly .. /u u- » » 1 oaa/i t » 1 ■^nAo^ d » 

,. ,r,^. , , tions (Hashimoto et al. 2004; Lamer et al. 2008). But, this ap- 

proportional to the baryon density squared. This makes clusters ■ , , j » n • 1 » 1 

^, . . ^ ■ , \r , proach has never been used systematically in cluster surveys, 

of galaxies prominent sources in the X-ray sky. • ..u » \. e • 1 * 

° ^ ^ ^ The main reason is that most of the existing cluster surveys 
Thanks to the spatial resolution of Chandra and the high sen- 
sitivity of XMM, coupled with their excellent spectral resolu- 1 „, ] T~ ,,. , c .1. ^ a t- t u u i j 
, , ., , , We note that a blind survey or the 6.4 Fe line has been applied 
tions, it has become possible to study the detailed thermodynam- X-ray surveys with interesting results (Maccacaro et al. 2004; 
ical properties of the ICM. In particular, the detection of emis- graito et al. 2005). The search of line emission in deep AGN surveys, 
sion lines of highly ionized metals, in both the core and outer however, is significantly different with respect to the search of Fe line 
emission in the thermal spectra of ICM, and their method is unsuitable 
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are based on source samples selected by ROSAT, whose energy 
range is 0.1-2.4 keV, hence does not cover the hard band where 
the Fe lines lie. To date, most ROSAT clusters have been con- 
firmed through optical imaging and spectroscopic observations. 
Several studies based on Chandra and XMM data (Tozzi et al. 
2003; Ettorietal. 2004; Balestra et al. 2007; Maughanetal. 
2008) simply confirmed the excellent agreement of X-ray and 
optical redshifts. This agreement was achieved using the optical 
redshift as the initial redshift parameter in the X-ray spectral fit. 
The blind search of the redshift from the X-ray spectral analysis 
has never been explored thoroughly. 

The situation could change significantly with the next- 
generation X-ray surveys of planned future missions sensi- 
tive in the 0.5-10 keV band, such as eROSITA (Predehl et al. 
2010) or the proposed Wide Field X-ray Telescope (WFXT, see 
Murray & WFXT Team 2010). In these surveys, we expect to 
detect several hundred thousands of groups and clusters of galax- 
ies, making an ad hoc follow-up program infeasible, and requir- 
ing careful coordination with existing and future optical and IR 
surveys (Rosati et al. 2010). Therefore, the ability to recoverred- 
shifts on the basis of the X-ray data alone, would be critical 
since it would allow one to measure the redshift, thus the in- 
trinsic physical properties of ICM, such as the X-ray luminosity 
Lx, ICM temperature Tx, iron abundance of the ICM Z, and 
the electron density Tig. In this case, it will be not only possible 
to investigate the physics of the ICM, but also perform cosmo- 
logical tests with large cluster surveys entirely characterized on 
the basis of X-ray data. 

The factors which affect the measurement of redshift from 
X-ray spectra are; the signal-to-noise ratio, hereafter S/N (or the 
total number of detected counts), the energy dependence of the 
instrument effective area, its spectral resolution, the intrinsic Fe 
abundance (or metallicity in general), the ICM temperature, and 
the actual redshift. A possible departure from collisional equilib- 
rium, or the incompleteness of the atomic models used in the fit- 
ting procedure, are not investigated here, and are expected to be 
negligible in the analysis of high-z cluster spectra. In this work, 
we use the data archive of the Chandra X-ray satellite to investi- 
gate the ability of recovering the X-ray redshift zx with a blind 
search of X-ray emission lines. This can lead to interesting ap- 
plications for future X-ray survey missions. 

The paper is organized as follows: in Section 2, we describe 
the data we use for our study. In Section 3, we present a sim- 
ple and efficient strategy to search for emission lines in X-ray 
spectra. In Section 4, we explore the use of weak priors to test 
their effect on X-ray redshift measurements. In Section 5, we 
apply our algorithm to the expected results from the WFXT 
surveys. Finally in Section 6, we summarize our conclusions. 
Throughout the paper, we assume the 7-year WMAP cosmology 
flA = 0.728, n,n = 0.272, and Hq = 70.4 km s'^ Mpc'^ 
(Komatsu et al. 2010). 

2. The data 

We searched the Chandra archive for targeted observations of 
galaxy clusters with different temperatures and redshifts. We 
adopted a sample of 46 clusters (see Table 1) mostly based 
on those studied by Balestra et al. (2007). This sample includes 
clusters with temperatures from ^ 4 keV to ^ 14 keV over a 
wide redshift range (0.16 < z < 1.4). The sample is incomplete 
and represents a collection of targets from different surveys (see 
reference list in Table 1). We also remark that in several cases 
we did not make use of all the available archival observations, in 
order to explore a wider S/N range. However, the most impor- 
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Fig. 1 : Detected photons as a function of redshift for the 46 clus- 
ters selected from the Chandra archive. The wide range of values 
reflects the different fluxes of each object and the different expo- 
sure times of the available data. 



tant point is that this sample covers a sizable portion of redshift- 
temperature space allowing us to comprehensibly explore the 
ability to measure the redshifts from the X-ray spectra. 

The number of detected photons in the 0.5-7 keV band 
ranges from ~ 230 to ~ 45000 as shown in Figure 1. Most 
were gathered with Chandra ACIS-I in the FAINT/VERYFAINT 
mode, while only six clusters were observed with ACIS-S. 
Calibration files were obtained using the most recent release of 
CALDB at the time of writing (CALDB 4.2). Image reduction 
began from the level 1 event file. We applied a charge trans- 
fer inefficiency (CTI) correction, FAINTA'ERYFAINT clean- 
ing, grade correction, and time-dependent gain correction. High 
background intervals were filtered with a 3a clipping procedure. 
The response matrix files (RMF) and the ancillary response files 
(ARF), necessary for the X-ray spectral analysis, were generated 
with mkacisrmf and mkwarf . For details about data reduc- 
tion, we refer to Balestra et al. (2007). 

The spectra were analyzed with XSPEC vl2.6.0 (Arnaud 
1996). Data were fitted with a single-temperature mekal model 
(Mewe et al. 1985, 1986; Kaastra 1992; Liedahl et al. 1995) in 
which the ratio of the elements was fixed to the solar value as in 
Anders & Grevesse (1989). Galactic absorption is modeled with 
tbabs (Wilms et al. 2000). The optical redshifts were collected 
from the literature. The typical error in the optical redshift is 
Az ^ 0.003. We adopt the optical redshift when measuring the 
temperatures and metallicities of the ICM. Our results, shown 
in Table 1, differ slightly from previous analyses presented in 
Balestra et al. (2007) because of the updated Chandra calibra- 
tion and in some cases the different extraction region. Computed 
errors always correspond to the Icr confidence level. 

3. Blind search for the redshift with no priors 

3.1. Method 

The Ka shell complex of Fe consists of two groups of lines. 
The first is the He-like iron (Fe XXV) K-shell complex, whose 
resonant line is at 6.7 keV, that is the most prominent spectral 
feature of ICM spectra at high temperatures (above 2 keV). The 
second group corresponds to the H-like iron (Fe XXVI) spectral 
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line whose resonant line is at 6.9 keV. Both groups include lines 
corresponding to several transitions as predicted in the case of 
collisional ionization equilibrium. The H-like and He-like com- 
plexes, separated by about 0.2 keV, can be resolved at the spec- 
tral resolution of Chandra ACIS in data of high S/N; in this case, 
the line ratio can be used as a temperature diagnostic in addi- 
tion to the continuum shape. On the other hand, the structure of 
each group is not resolved and reflects the asymmetric shape of 
the observed feature. Both effects are taken into account by the 
mekal model, which includes excitation, radiative recombina- 
tion, dielectronic recombination, and inner shell ionization, and 
assumes the optically thin limit (i.e., that no photo-ionization or 
photo-excitation effects are taken into account). Given the low 
S/N of distant cluster spectra, the line complex often appears 
as a single, prominent feature. 

Owing to uncertainties in the ACIS caUbration below 0.5 
keV, we excluded these low energy photons from the spectral 
analysis to avoid any systematic bias. The effective cut at high 
energies is generally around 7 keV, since the S/N of a thermal 
spectrum rapidly decreases above 5 keV. Therefore, we consider 
the energy range 0.5-7 keV for both imaging and spectral analy- 
sis. Given the explored redshift interval (0 < z < 2), the Fe line 
complex is always well within this range. 

Background is usually selected in the same chip where 
the source lies. When several observations are used, the back- 
ground can only be selected in overlapping regions. For some 
nearby clusters that occupy a whole chip, or even an entire 
field (Abell 1689 for example), we compile a synthetic back- 
ground using CALDB. The centroid of the cluster emission 
is determined by surface brightness fits adopting a standard 
/3 model (Cavaliere & Fusco-Femiano 1976). A set of circu- 
lar regions are then drawn to select the extraction radius Rext 
that maximizes the S/N, computed as SN{R) = {S{R) - 

B)/^S{R) + {As/Ab) X B, where S{R) are the total num- 
ber of photons detected within the radius R, and B are the back- 
ground number of photons expected in the same area. The factor 
As/Ab is the geometrical backscale, i.e. the ratio of the area of 
the source extraction region to that of the background. Since our 
sources are extended and the background is chosen to be as close 
as possible to the source, this factor can be of order < 1. This ex- 
pression properly takes into account the statistical uncertainty in 
the number of photons in the source region and the region where 
the background is sampled. The extraction regions are selected 
on the X-ray image in the total band (0.5-7 keV). 

We adopt the default Levenberg-Marquardt fitting algorithm. 
There are four parameters to be fitted: temperature, metal abun- 
dance, redshift, and normalization. To test the capability of re- 
covering the actual redshift with a blind search of the Kq, Fe line 
in the X-ray spectrum, we repeat our fits starting from a refer- 
ence set of parameters, corresponding to Zstart ~ 0.1, kTstart = 
5 keV, and Zstart /Zq — 0.3 in units of Anders & Grevesse 
(1989) for all the clusters. We use Cash statistics (Cash 1979) ap- 
plied to an unbinned source plus background counts, and there- 
fore exploit the full spectral resolution of the ACIS -I and ACIS- 
S instrument. Cash statistics ensure better performance with re- 
spect to the canonical analysis of binned data, particularly for 
low S/N spectra (Nousek & Shue 1989). After finding the abso- 
lute minimum, we explore the redshift space with the steppar 
command, covering the entire range of possible values from 
z = 0.01 to z = 2.0 with a very small step Sz = 0.01. When a 
new minimum is obtained, the best fit is automatically updated. 
We then plot the difference of the C-stat value with respect to the 
minimum as a function of redshift. The ACstat{z) function is 
therefore the difference between the absolute minimum and the 



best-fit value obtained when the z parameter does not match any 
line, which is equivalent to optimizing the fit when the Fe abun- 
dance is forced to be zero. One example (MACSJ1423) is shown 
in Figure 2. The ACstat rapidly declines to a minimum when- 
ever a Fe line candidate is found. The deepest minimum in the 
C-statistics as a function of the redshift in Figure 2 clearly shows 
that the best-fit redshift agrees with the optical redshift (indi- 
cated by the vertical line). The horizontal line, therefore, cor- 
responds to the minimum Cstat value obtained with zero metal 
abundance. 

The case of a catastrophic error is shown in Figure 3. Here 
there are no minima corresponding to the optical redshift. We 
note that, instead, there are several comparable secondary min- 
ima. We explore a possible use of secondary minima in Section 
3.7. 



Fig. 2: ACstat as a function of the redshift parameter for 
MACSJ1423. The vertical line indicates the optical redshift. The 
horizontal line corresponds to the minimum Cstat allowed when 
the metal abundance is set to zero. 



3.2. Results 

For most clusters, the redshifts zx found from the X-ray spec- 
tral analysis are consistent with their optical values Zo within 
the errors as shown in Figure 4. However, there are still sev- 
eral clusters for which zx is far from Zq. To define these "catas- 
trophic failures", we perform a Sa-clipping by computing the 
rms value Az^ms = ■\/S(zx — Zo)^/n for the entire sample of 
46 clusters. The catastrophic errors are shown as empty circles in 
Figure 4 and as an empty histogram in Figure 5, where we plot 
the quantity Az / Az„ns- Thanks to this definition, we identify 
eight catastrophic errors. After removing the catastrophic fail- 
ures, the rms redshift deviation is Az ~ 0.03. We note that this 
uncertainty is somewhat larger than the typical statistical error 
az estimated by the XSPEC fitting routine. We argue that the 
statistical error in the redshift is slightly underestimated. As we 
see later (Section 5.1), ACstat is not only the most relevant pa- 
rameter to accept or reject the zx value, but can also be used to 
adjust the statistical eiTor obtained by the spectral fit thanks to 
its strong correlation with the actual value of zx — Zq. 

Our results for the entire sample of Chandra clusters is 
shown in Table 2. From a comparison with Table 1, it is pos- 
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Fig. 3: ACstat as a function of the redshift parameter for 
SCI 120. The vertical line indicates the optical redshift. At vari- 
ance with Figure 2, the search for a minimum in Cstat leads to a 
catastrophic error. 



Fig. 5: Histogram of the quantity Az/Az^ms- Catastrophic fail- 
ures, identified with a 3ct clipping, are shown with an empty 
histogram. A Gaussian fit to the data is shown as a dashed line. 
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Fig. 4: X-ray redshift zx compared to the optical values Zq for 
the 46 Chandra clusters in our sample. The eight empty circles 
are the catastrophic failures selected with a 3a clipping. 



sible to see that, apart from the catastrophic failures, the best-fit 
temperatures are unaffected, while the best-fit metal abundances 
are systematically higher (see discussion in Section 3.5). 

3.3. Tuning the energy band 

The Ka Fe Une is always located in the 2-7 keV energy range 
(observing frame). Therefore, before analyzing in detail our re- 
sults, we investigate whether focusing on this band may only 
help in driving the best-fit redshift value more efficiently towards 
the real one. We can either select the extraction region from the 
hard band image, or fit the spectrum in the hard energy band 
only, or adopt both criteria. The hard-band selected regions are 
usually smaller since the bulk of the X-ray photons from the ICM 
thermal spectrum are in the soft (0.5-2 keV) band. To compare 
the results of different strategies, we introduce a quality param- 



Region 


Energy band 
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^Zrrn-s 


# of cat. errors 


Total 


Total 


133.60 


0.0288 
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Total 


Hard 


252.79 


0.0378 
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Hard 


Total 


284.14 


0.1736 
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Hard 


Hard 


468.66 


0.0356 
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Table 3: Quality parameter for different choices of the energy 
band and typical rms value of Az after applying the 3a clipping. 



eter defined as Q = J^ii^x.i — Zo.i)l<^zxY' 1"^^ where the sum is 
performed for the entire sample (46 clusters) without removing 
the catastrophic failures. The value of Q is a useful estimate of 
the average discrepancy between zx and Zo obtained with differ- 
ent algorithms. The comparison is shown in Table 3, where we 
also list the Az^ms values after applying 3a clipping, and the 
number of catastrophic failures. We find that our default choice 
based on the total energy band (0.5-7 keV) has a clear advantage 
with respect to strategies focusing on the hard band. We argue 
that the signal in the soft band is useful for defining at best the 
continuum, and this, in turn, positively affects the detection of 
the Fe line. We therefore conclude that there is no gain in re- 
stricting the energy band to the 2-7 keV range, and that the full 
amount of spectral information is useful in searching for the Fe 
Kq line. 



3.4. Dependence on the net detected photons 

We investigate how the average deviation and the number of 
catastrophic failures depend on the S/N of the spectra or, alter- 
natively, on the total number of net detected photons. In Figure 
6, we show the quantity zx — Zo as a function of the net de- 
tected photons. From the inspection of Figure 6, we note that it 
is unclear whether a threshold on the number of detected photons 
exists above which the X-ray redshift can be considered reliable 
within a given confidence level. When we consider S/N instead, 
a similar result is obtained. We note that above a threshold of 
about 1000 counts (vertical line) the number of catastrophic fail- 
ures is not negligible, although five out of eight catastrophic fail- 
ures are below this threshold. Three catastrophic errors occur in 
spectra with thousands of counts. We draw the conclusion that 
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Fig. 6: Difference between optical and X-ray redshifts as a func- 
tion of the detected photons. The eight empty circles are catas- 
trophic failures. In some cases, spectra of high S/N (e.g. with 
more than 1000 net photons, shown with a vertical line) may not 
be successfully fitted. 



Fig. 7: "True" Fe abundance (obtained by fitting the X-ray spec- 
tra freezing the redshift to the optical value) versus the total num- 
ber of detected photons. The radius of the dots is proportional to 
I Az|. The three empty circles represent 1 a upper limits. 



the efficiency of measuring the redshift depends not only on the 
high quality of the signal, but also on the intrinsic properties of 
the source. The likelihood of obtaining a reliable spectral char- 
acterization, including the redshift, does not depend only on the 
net number of photons. To find a robust method to select reliable 
X-ray redshifts, we proceed with a deeper investigation of the 
spectral parameters. The most obvious parameter to consider is 
the actual Fe abundance. 

3.5. Dependence on the intrinsic metal abundance 

We explore now whether the intrinsically low Fe abundance can 
be another relevant source of error. In Figure 7, we plot the "true" 
Fe abundance Zpe (i-S-, the value obtained when the redshift 
is fixed to the optical value, see Table 1) versus the number of 
net detected photons in each spectrum. Most catastrophic fail- 
ures are located in the bottom-left corner. This indicates that the 
largest discrepancies are associated with both the low S/N of 
the X-ray spectrum and the low Fe abundance. Significant dis- 
crepancies can be found among spectra with more than ^ 4000 
photons if the Fe abundance is Zpe < 0.2Zq. 

This information cannot be used in our blind search for the 
emission lines, since we do not know a priori, the Fe abundance. 
As we see in the next section, making an assumption about the 
intrinsic Zpe value does not improve the fit. In any case, we 
find that it is impossible to use a sample whose redshift is de- 
termined via the X-ray analysis to investigate the Fe abundance. 
The X-ray fits with a free redshift parameter will always lead 
to values that are systematically higher than those obtained by 
fixing the redshift to the optical value Zq, as shown in Figure 8. 
We find that the typical value of the ratio of the fitted Fe abun- 
dance to the "true" one is {Zpe{zx)/Zpe{zo)) = 1.110, which 
corresponds to a positive bias of 11%. This result is expected 
since the position of the Fe line is found by maximizing the Fe 
abundance for a given temperature. On the other hand, there is 
no evidence of a bias in the best-fit values of the temperature, 
as shown in Figure 9. For the temperature ratio, we find that 
{T{zx)/T{zo)) - 0.993. 




Fig. 8: "True" Fe abundances Zpeizo) in the Chandra sample 
(obtained by fixing the redshift parameter to the optical value) 
plotted against those measured leaving the redshift parame- 
ter free. Catastrophic failures are not included. The values of 
Zpeizx) are systematically higher than those we derive when 
the redshift is fixed to the optical value. 



3.6. Dependence on ACstat 

Another important indication is the intensity of the emission 
line itself. The most reliable lines are those that provide the 
largest decrease ACstat in the C-statistics (see Figure 2). A 
larger ACstat indicates a more robust emission line, then a more 
reliable X-ray redshift, as shown in Figure 10. Assuming that 
ACstat behaves similarly to the Ax^ for one degree of freedom, 
we find that a threshold that excludes all the catastrophic errors 
is ACstat > 9, corresponding to a nominal confidence level of 
3(7. Above this level, we have 26 clusters for which the redshift 
is measured with good accuracy, constituting a sample virtually 
free of catastrophic errors. 
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Fig. 9: "True" temperature T{zo) (obtained when the redshift 
parameter is fixed to the optical value Zo) plotted against the 
value T{zx) measured leaving the redshift parameter free. 
Catastrophic failures are not included. 
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In summary, we do not find any benefit when considering 
the secondary minima, and we confirm that simply using a sharp 
threshold on ACstat is the most efficient criterion for selecting 
reliable redshift measurements. 



4. Refining thie strategy with weal< priors 

We have previously here applied a direct blind search for the Fe 
line without making any further assumption. We now investigate 
whether we can achieve a more efficient strategy using weak pri- 
ors. Two possible choices of priors come from the study of local 
cluster samples: the low scatter in the measured Fe abundance of 
hot clusters and the tightness of the L — T relation. 

4. 1. Prior on the Fe abundance 

The Fe abundance is observed to be almost constant in local and 
medium redshift clusters for virial temperatures kT > 5 keV 
(Renzini 1997; Baumgartneret al. 2005; Balestra et al. 2007). 
This suggests that we may remove one fitting parameter by 
freezing the Fe abundance to Zpe = 0.3Zq. This assumption 
would be wrong below 3 ke V, since in the low temperature range 
(from poor clusters to groups) the Fe abundance values cover 
a wide range (Renzini 1997; Rasmussen & Ponman 2007). To 
explore the effects of this assumption, we repeat the fits with 
Z = 0.3Zq and show the results in Figure 11. The Q value 
after this assumption is 103.12, which represents a mild im- 
provement in accuracy, mainly because of the larger error bars of 
MACSJ0404 and RXJ0542. Incidentally, we note that these two 
clusters with "tnie" Fe abundance below 0.2Zq are fitted with a 
redshift ^ 2.1, a mistake caused by an edge in the Chandra re- 
sponse around 2 keV. The average rms value of Az is also 0.03, 
and the number of catastrophic failures is 9. We conclude that 
there is no evidence that a prior on the Fe abundance provide 
more reliable results. In addition, we find some indication that 
this assumption might prevent us from being able to use the cri- 
terion about ACstat to select the most reliable redshift measure- 
ments, as shown in Figure 12, where we have two catastrophic 
failures for ACstat > 9. 



Fig. 10: The quantity Az = zx — Zo plotted against ACstat - The 
eight empty circles are catastrophic failures. The vertical line 
show the nominal 3ct threshold corresponding to ACstat = 9. 



3. 7. Exploring tlie secondary minima 

We have investigated the properties of the absolute minima. In 
principle, the Cstat function of each cluster has several sec- 
ondary minima that may contain useful information. We ex- 
plore this possibility by identifying all the local minima in the 
Cstat — z function, and considering the corresponding ACstat 
and Az. We find that, as for the catastrophic failure, none of their 
secondary minima provide a correct redshift measurement. As 
for the 38 clusters with acceptable redshift measurements, sec- 
ondary minima do not provide any improvement in the redshift 
value, as expected. 

We also explore the possibility of flagging unreliable zx 
from the presence of one or more secondary minima close to the 
absolute minimum. However, applying a naive rejection based 
on more than one line of similar ACstat, would result in the re- 
jection of several successful cases. 



^ RXJ0542 




0.0 0.5 1.0 1.5 2.0 

Fig. 11: Same as Figure 4, after assuming a fixed Fe abundance 

Zpe = O.SZq. 
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Fig. 12: Same as Figure 10, after assuming a fixed Fe abundance 

Zpe = O.3Z0. 

4.2. Priors on the L - T relation 

Additional information is the clearly defined L — T relation 
measured for local groups and clusters of galaxies (see, e.g., 
Lumb et al. 2004; Maughan et al. 2006; Pratt et al. 2009). Large 
errors in the redshift will result in a position in the L ~ T plane 
significantly distant from the average observed relation, once the 
intrinsic scatter is properly taken into account. 

To illustrate this effect, we compute the bolometric luminosi- 
ties for our sample of clusters. The luminosity contribution out- 
side the extraction regions is estimated by fitting a (3 model to 
the observed surface brightness. All luminosities are computed 
for i?5oo, which is defined as the radius within which the average 
density contrast is 500 times the critical density of the universe 
at the cluster redshift. As for the L-T relation, we adopt the best 
fit by Maughan et al. (2006) to the Wide Angle ROSAT Pointed 
Survey (WARPS) sample at 0.6 < z < 1.0: 

isoo = 4.97 ± 0.80 i-^j W Vo' ergs"! . (1) 

In Figure 13, we compare this relation with that obtained 
with our cluster sample in a similar redshift range. Catastrophic 
failures are labeled with empty circles. In total 14 clusters, in- 
cluding six catastrophic errors, are formally excluded if we re- 
quire all the clusters to be consistent at the 3(7 level with the as- 
sumed L-T relation. We note, indeed, that the catastrophic errors 
are generally distant from the measured L-T relation, but still 
very close to the bulk of the cluster, so that a clear separation is 
not observed. This happens because the fitted temperature scales 
as (1 + zx), and the relation is approximately L oc T^. This 
implies that a wrong redshift would move the cluster luminosity 
approximately along the observed L ~ T relation. Another less 
relevant constraint that we show in Figure 1 3 is an upper limit 
on the ICM temperature, which is conservatively put to 20 keV 
(corresponding to the 3 — 8 x IQ^^Mq range). This constraint is 
expected to be useful in only a few extreme cases. 

In principle, this prior should be applied at different redshifts 
because the L — T relation and its intrinsic scatter are expected 
to evolve at some level. However, there is still significant un- 
certainty in the measured L-T evolution (see Ettori et al. 2004; 



Fig. 13: L-T relation for our cluster sample compared with the 
relation measured by Maughan et al. (2006). Catastrophic fail- 
ures are marked as empty circles. The red dashed lines indicate 
the 3(7 confidence levels around the average L-T relation. The 
vertical dashed line marks a reasonable upper limit (20 ke V) that 
can be assumed for ICM temperatures. 

Maughan et al. 2006; Branchesi et al. 2007). To summarize, the 
L-T prior can be effective given that slope, normalization, and 
intrinsic scatter are accurately known at different redshifts. In 
practice, this kind of information is expected to be obtained from 
cluster samples as large as those achievable with future wide 
area surveys. In this case, this criterion would not be a prior, but 
rather a self-calibration procedure. Therefore, with the present 
knowledge this prior should be applied with caution. 

5. Building a compiete sampie of dusters with 
measured zx - appiication to future X-ray 
surveys 

The work presented so far shows the efficiency of a blind search 
for a redshift by means of the X-ray spectral analysis of known 
clusters observed by Chandra. On the basis of these results, we 
can define general criteria for measuring accurate redshifts when 
applying this procedure to a sample of new clusters, whose red- 
shift is not known. We can apply this to new extended sources 
detected in the Chandra and XMM archives. However, serendip- 
itous searches in Chandra and XMM pointings, which are in- 
evitably confined to relatively small sky areas, result in detec- 
tions with a Umited amount of net counts, such as the ChaMP 
Galaxy Cluster Survey (Barkhouse et al. 2006), the XMM- 
Newton Distant Cluster Project (XDCP, Fassbender 2008), the 
XMM-Newton Large-Scale Structure (XMMLSS, Pierre et al. 
2007), and the XMM Clusters Survey (XCS, Romer et al. 2001; 
Lloyd-Davies et al. 2010). In this case, the application of our 
procedure will still be helpful but have a modest impact. 

The situation radically changes if we consider future X- 
ray surveys. The planned eROSITA satellite (Predehl et al. 
2010; Cappelluti et al. 2010) and the proposed WFXT mis- 
sion (Murray & WFXT Team 2010) will provide a large num- 
ber of new detections. We focus here on specific expectations 
for WFXT, which is the most optimized mission for surveys 
thanks to its wide-field optical design. This mission will be able 
to provide a golden sample of 15000-20000 rich clusters, with 



7 



H.Yu et al.: Measuring redshifts using X-ray spectroscopy of galaxy clusters 



kT > 3 keV out to z ~ 1.5, detected with more than 1500 
counts (Borgani et al. 2010). 

5. 1. Selecting clusters with reliable zx 

We use the simulations we performed to assess WFXT science 
cases as described in Tozzi et al. (2010). Using the halo mass 
function by Sheth & Tormen (1999), we first determine the cos- 
mological parameters that most accurately fit the observed X- 
ray cluster luminosity function (XLF, see Rosati et al. 1998; 
Giacconi et al. 2002; Bauer et al. 2002) for a given relation be- 
tween cluster mass and X-ray luminosity. In this way, we can 
extrapolate the observed XLF over the luminosity and redshift 
ranges expected to be covered by the WFXT surveys. This lumi- 
nosity function is then used to generate mock realizations of the 
WFXT cluster surveys through a Monte-Carlo procedure. The 
corresponding virial temperature (temperature gradients are not 
considered) is assigned with a scatter of 10% with respect to 
the mass, and the luminosity with a scatter of 15% - 45% (going 
from massive clusters to groups) with respect to the temperature, 
to be consistent with the observed L-T relation of local clusters. 
We are thus able to simulate an accurate spectrum normalized 
to the predicted flux of each cluster Metal abundances are as- 
signed randomly with a Gaussian distribution centered around 
the average value Z = 0.3Zq with a sigma AZ = O.IZq, with 
a hard lower limit at Z = O.IZq. Values below Z — O.IZq 
have never been measured in hot clusters at any redshift. In ad- 
dition, stacked spectra analyses of distant clusters identify sig- 
nificant enrichment in the ICM up to z ~ 1.3 (Balestra et al. 
2007; Maughan et al. 2008). As a matter of fact, a safe assump- 
tion would be a constant Fe abundance equal to 0.3^0, as also 
suggested by Figure 7. However, we believe that by allowing 
Zpe values as low as O.IZq, we will provide a conservative es- 
timate of the number of line detections. We also note that with 
this choice, in the simulation analysis we obtain several best-fit 
values of the Zpe parameter populating the range below O.IZq. 
At present, it is impossible to provide a more accurate model for 
the distribution of Zpe in the ICM as a function of redshift. 

The background is computed accordingly to expectations for 
the WFXT mission. The net count rates per deg^ of the different 
background components are given in Table 4. We assume that the 
source spectrum is extracted from a 500 kpc region, a value cal- 
ibrated on the analysis of real Chandra data as in Balestra et al. 
(2007). We also take into account that typically 10% of the total 
emitted flux is lost outside the extraction regions. The effect of 
vignetting is also considered, since the clusters will be randomly 
positioned across the field of view. This is an important factor 
to be included because the vignetting is particularly severe in 
the hard band, where the Fe line is located. We use effective area 
files corresponding to seven different off-axis angles in the range 
0-30' to reproduce as closely as possible the vignetting effect on 
the observed spectra. 

We produce a mock catalog of groups (with temperatures 
above 0.5 keV) and clusters as expected in 100 square degrees 
of the WFXT Medium Survey (see Rosati et al. 2010). The sim- 
ulation thus consists of 100 pointings of 13.2 ks each, extracted 
from a total 3000 square degree area. We have about 2500 groups 
and clusters above the flux of 4 x 10^^^ erg s^^ cm^. For all of 
them, we measure the redshift zx with our procedure described 
in Section 3. We remark that this mock sample, being flux lim- 
ited, includes a large number of clusters with kT < 2 keV, for 
which the presence of the L-shell line complex at low energy 
significantly facilitates the measurement of zx- 



Background component 


0.5-2 keV 


2-7 keV 


Particles 


0.188 


0.397 


Galactic 


21.4 


0.0 


AGN (Medium survey) 


3.9 


1.65 


Cluster (Medium survey) 


0.79 


0.14 


Total (Medium survey) 


26.28 


2.19 



Table 4: WFXT background net count rates per deg^ in the 
Medium Survey (13.2ks exposure). The AGN and cluster contri- 
butions correspond to the emission of the unresolved AGN and 
groups/clusters, respectively, in the Medium Survey. 
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Fig. 14: Percentage of catastrophic errors as a function of 
ACstat in a sample with different net photons within Rg from 
top to bottom. The horizontal line marks the 5% level. 



The catastrophic errors are identified by means of 3cr clip- 
ping. In Figure 14, we plot the percentage of catastrophic fail- 
ures as a function of ACstat for clusters detected above a given 
photon threshold. We see that, for a given ACstat, the number 
of catastrophic failures rapidly increases below 1000 net pho- 
tons. By selecting a reference sample with more than 1000 net 
photons and ACstat > 9, we are able to maintain the percent- 
age of catastrophic errors below 5%. The effect of the crite- 
rion ACstat > 9 can be appreciated by comparing Figs 15 and 
16. In Fig. 15, we show all the groups and clusters with more 
than 1000 net photons, while in Fig. 16 we apply the threshold 
ACstat > 9, which helps to define a "golden" sample of 862 
groups and clusters. This sample still includes 38 catastrophic 
errors (4.4% of the total) but all of them are rejected marginally 
with redshift deviations slightly larger than 0.032. After apply- 
ing the cut for ACstat > 9, the mean redshift offset (Az) is 
-0.0024 and the deviation rms Azrms is 0.01 16 (see Figure 17). 
We note that, despite 38 measurements being classified as catas- 
trophic errors, their typical discrepancy is still very small, and 
the overall redshift accuracy of the sample with ACstat > 9 sat- 
isfies the demands of cosmological tests. For example, in order 
not to degrade dark energy constraints by more than 10%, both 
(Az) < 0.003 and Az^ms < 0.03 are required (Lima & Hu 
2007), which are satisfied by our "golden" sample. 

Figure 18 shows Az ~ zx — Zq as a function of ACstat for 
all the clusters with detections of more than 1000 photons. The 
error rapidly decreases with increasing ACstat, as expected. The 
relation between Az^ms and ACstat can be approximated as 
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Fig. 15: X-ray measured redshift zx for all the clusters in a 100 
square degree field of the WFXT Medium Survey detected with 
more than 1000 photons within the extraction region. We have 
1037 clusters and groups, and 118 catastrophic failures. The col- 
ors are set according to the cluster temperature. 



Fig. 17: Histogram of Az/Az^ms from a simulated 100 square 
degrees of the WFXT Medium Survey. The sample here include 
all clusters with more than 1000 photons and ACstat larger than 
9. Empty histogram indicates the residual catastrophic eiTors. 
The red dashed line shows the Gaussian fit after 3a clipping. 




1.0 

True Redshift 




Fig. 16: The same as Figure 15, but after applying the criterion 
ACstat > 9. There are no deviations larger than 0. 1 , and only 38 
clusters with Az > 0.032. The size of the sample has decreased 
to 862. 



Fig. 18: Az ~ zx^Zq versus ACstat for clusters with more than 
1000 photons. Error bars show la statistical error from XSPEC 
fitting. The color code is the same as in Figure 15. There is an 
obvious trend of Az decreasing at larger ACstat- The red solid 
lines delineate la envelope for Az„„s- The red dashed lines de- 
lineate the 3(7 confidence levels. 



Az„„, = OmbilogwACstat)-^ '' ■ (2) 

With this information we can select a sample with reliable red- 
shifts and a robust error estimate for zx- 

5.2. Assessing the completeness of the sample 

Samples selected with source counts above a given number of 
net photons can still be treated as "complete". Such a thresh- 
old corresponds to a good accuracy for the limiting flux as a 
function of the position on the sky, once the exposure time and 
the vignetting in each field of view is computed. Therefore, it is 
possible to obtain a clearly defined sky coverage from which the 



XLF and the X-ray Temperature Function (XTF) can be com- 
puted for cosmological applications. 

However, it is clear that a criterion based solely on the net 
counts is an insufficient one. As is immediately visible from 
Figure 15, among the clusters whose X-ray spectra have more 
than 1000 photons, there are about 120 catastrophic failures, 
several of which are at a redshift much higher than the true 
one. This implies that there may be a significant contamination 
at the bright ends of the luminosity and temperature functions, 
which are both particularly sensitive to cosmological parame- 
ters. For example, six clusters in our sample with more than 
1500 counts are incorrectly located at z > 1.5. For the full 
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3000 deg^ Medium Survey, this would correspond to about 180 
fake z > 1.5 clusters, compared to the 20 expected (see Figure 
2 in Borganiet al. 2010). This implies that the majority of the 
rarest z > 1.5 clusters candidates would be spurious, and there- 
fore that the constraints on both cosmological parameters and 
non-Gaussianity (e.g. Sartoris et al. 2010; Verde 2010, and ref- 
erences therein) would be highly biased without applying the 
^Cstat > 9 criterion. 

By adopting the threshold on ACstat, the price to pay is 
a lower level of completeness of the sample. In the aforemen- 
tioned simulation, we have excluded 175 clusters (17% of the 
total) detected with a number of photons above the threshold 
because they do not satisfy the ACstat > 9 criterion. In this 
way, we remove most of the catastrophic errors but also lose 
about 100 clusters with a good z^. As a result, the "golden sam- 
ple" is incomplete above a given flux threshold, and the missing 
clusters must be accounted before applying cosmological tests. 
Unfortunately, this step is difficult because the temperature dis- 
tribution of the clusters with ACstat < 9 is biased with respect 
to the distribution of the entire sample. This is shown in Figure 
19, where the normalized distribution of rejected clusters with 
more than 1000 photons is skewed towards high temperatures. 
This occurs because it is generally more difficult to fit the Fe 
line in increasingly hot clusters, since the abundance of He-like 
Fe ions, which dominate the line emission, begins to decrease at 
temperatures kT > 5 keV, and the thermal continuum steadily 
increases. Such a bias against the most massive objects needs to 
be carefully quantified. 

There are two ways to restore the completeness of the sam- 
ple. One is to statistically correct the incompleteness using re- 
sults from N-body simulations. However, this is affected by se- 
rious systematic errors because we modify the high-mass end 
of the cluster distribution, which is the most sensitive to cos- 
mological parameters. The other option is to limit spectroscopic 
foUow-up to only those clusters with ACstat < 9. 

To design the optimal strategy, we should proceed with a 
detailed investigation of the effects of varying the threshold of 
ACstat on the cosmological parameter constraints. We noted 
in Section 5.1 that ACstat > 9 provide a "golden sample" 
that largely satisfies the requirements of achieving a precision 
of 10% for the dark energy constraints. Relaxing the threshold 
on ACstat and mapping the threshold values to the accuracy on 
the cosmological parameter is clearly an option that we plan to 
explore in a future paper. 

Another possible source of contamination may be diffuse 
X-ray emissions produced by inverse Compton processes as- 
sociated with radio jets (see, for example, Fabian et al. 2003). 
These sources are a serious contaminant when detecting X- 
ray extended sources, especially at very low fluxes. However, 
this source will hardly provide a spurious line detection above 
the selection threshold, hence their presence in a cosmological 
"golden" sample can be neglected. 

6. Conclusions 

We have performed a blind search of the Fe line in the X-ray 
spectra of clusters of galaxies using Chandra archival data. Our 
goal is to define the optimal methodology for determining X-ray 
redshifts and to quantify their accuracy by investigating both sta- 
tistical and systematic errors. To this end, we compared zx with 
the value Zo determined with optical spectroscopic observations. 

We have found good agreement in general, but also in sev- 
eral cases where zx is determined by the false detection of 




Input T(keV) 



Fig. 19: Temperature distribution for the sample of groups and 
clusters detected with more than 1000 net photons and ACstat > 
9 (filled histogram). The black line indicates the complete in- 
put catalog, while the red histogram refers to the clusters re- 
moved from the flux-limited sample after applying the criterion 
ACstat > 9. 



an emission line. Thus we explored different methods to min- 
imize the number of these catastrophic failures. We have found 
that it is preferable to use the total available band (0.5-7 keV) 
for both selecting the extraction region and fitting the spectra. 
Catastrophic failures are found to be caused by low spectrum 
S/N and an intrinsically low Fe abundance. For this reason, a 
lower limit to the net detected photons is insufficient to guaran- 
tee a robust measurement of the redshift. However, we find that 
we can exclude most catastrophic errors by using the criterion 
ACstat > 9, where ACstat is the difference between the best-fit 
value of Cstat, and its best-fit value when the metal abundance 
is constrained to be zero. 

We have also explored the use of weak priors to improve the 
results. We have found that by fixing the Fe abundance to the 
local average value Zpe ~ O.SZ©, the number of catastrophic 
errors is not reduced. Furthermore, by requiring that all the fitted 
clusters agree with the observed local L-T relation and its scat- 
ter, one has an independent method to identify outliers, although 
not as efficiently as the condition ACstat > 9. 

Our spectral analysis of Chandra clusters shows that future 
X-ray survey missions will be able to define sizable samples 
of clusters with X-ray measured redshifts. We specifically in- 
vestigated the case of the Wide Field X-ray Telescope using a 
mock catalog of ^ 2500 groups and clusters extracted from a 
100 square degree area of the WFXT Medium Survey. The input 
cluster catalog is generated from the cluster mass function nor- 
malized to the observed space density of X-ray clusters by as- 
suming local scaling relations to derive X-ray luminosities and 
temperatures associated with a given mass. We found that, by 
applying the condition ACstat > 9, we can successfully mea- 
sure the redshift of 862 clusters with more than 1000 net counts, 
which leaves a small fraction (4.4%) of marginally catastrophic 
errors. We argue that these subsamples with X-ray determined 
redshifts can be effectively used for cosmological applications, 
thus avoiding time-consuming spectroscopic observations, al- 
though additional simulations will have to be developed to assess 
the completeness of the sample down to a given flux limit. 
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ID 


Zo 


ref 


ObslD 


Detector/Mode 


texp (^'^) 


T{zo){keV) 


ZFe{Zo)/ZQ 


nHicm '') 


AbelI907 


0.160 


1 


3185 


1-V 


47 


5.56_Q QQ 




5.65x10^" 


AbelI1689 


0.187 


2 


540,1663 


1-F 


21 


13.74+°«^ 




1.82x10^° 


1E0657-56 


0.296 


3 


554 


1-F 


25 


is.nt°oil 


n 9+0.06 
—0 06 


6.5x10^" 


MS1008. 1-1224 


0.306 


4 


926 


1-V 


44 


6-09+0 33 


r, 97+0.06 
^•■^ ' —0.06 


7.26x10^" 


MS2137.3-2353 


0.313 


5 


928 


S-V 


34 


5.26«-i? 




3.55x10^° 


Abelll995 


0.319 


6 


906 


S-F 


56.4 


q 1O+0.32 

— 31 




1.42x10^" 


ZwC11358.1-f6245 


0.328 


8 


516 


S-F 


48.3 


7 -.9+0:26 
■ —0 25 


0.39+°°^ 


1.92x10^" 


MACSJ0404.6+1109 


0.355 


9 


3269 


1-V 


21.6 




0.16±°'?I 


1.43x10^° 


RXJ0027.6+2616 


0.367 


9 


3249 


1-V 


9.8 


8.07+}-g 


0.52+^,11 


3.86x10^° 


MACSJ1720.2+3536 


0.387 


10 


3280 


1-V 


20.8 






3.4x10^° 


ZwC10024.0+1652 


0.395 


11 


929 


S-F 


39.5 


4.57l°|? 


0.75i°'?^ 


4.23x10^° 


RXJ1416.4+4446 


0.400 


12 


541 


1-V 


31 


3.56_Q 




1.22x10^° 


MACSJ0159.8-0849 


0.405 


13 


3265 


1-V 


17.6 






2.08x10^° 


MACSJ2228.5+2036 


0.412 


9 


3285 


1-V 


20 


9r+0.61 


0.35iro9 


4.58x10^° 


MS0302.7+1658 


0.424 


8 


525 


1-V 


10 


4.38in^ 




10.9x10^" 


MACSJ0417.5-1154 


0.44 


14 


3270 


1-V 


12 




0.29tVA 


3.86x10^° 


MACSJ 1206.2-0847 


0.44 


15 


3277 


1-V 


23 






3.72x10^° 


RXJ1701. 3+6414 


0.453 


16 


547 


1-V 


49 


4.49l°-^g 




2.46x10^" 


RXJ 164 1.8+4001 


0.464 


16 


3575 


1-V 


45 


4.8llO-g 




1.1x10^° 


MACSJ 1824.3+4309 


0.483 


10 


3255 


1-V 


14.8 


SMtlil 


< 0.23 


4.58x10^° 


MACSJ13 11.0-0311 


0.494 


17 


3258 


1-V 


14.8 


8.671^1? 
5.58_Q 


0-37lVil 


1.87x10^" 


RXJ1524.6+0957 


0.516 


16 


1664 


1-V 


50 




2.91x10^° 


MS0015.9+1609 


0.541 


8 


520 


1-V 


67 


q+0.32 

— 32 


0.3liro5 


7.26x10^° 


MACS J 1423. 8+2404 


0.543 


18 


1657 


1-V 


18.5 




3+°'i 


2.38x10^° 


MACS J 1149.5+2223 


0.544 


18 


1656, 3589 


1-V 


38 


12.75_Q gr^ 


0.27+°'i 


2.28x10^° 


SCI 120-1202 


0.562 


7 


3235 


I-V 


68 


6.49+}?'^ 


0.21±°'f9 


5.19x10^° 


MS2053.7-0449 


0.583 


8 


551, 1667 


1-V 


88 






4.96x10^" 


RXJ0956.0+4107 


0.587 


16 


5294 


1-V 


17.2 






1.14x10^° 


MACSJ2129.4-0741 


0.589 


18 


3199 


1-V 


17.6 


9.48lJ-f 




4.86x10^° 


MACSJ0647.7+7015 


0.591 


18 


3196 


1-V 


19.2 




<0.1 


5.64x10^" 


RXJ0542.8-4100 


0.634 


7 


914 


1-F 


50 






3.73x10^° 


MACSJ0744.9+3927 


0.698 


18 


3197, 3585 


1-V 


40 




0.29t^-ll 


5.7x10^° 


RXJ 122 1.4+49 18 


0.70 


12 


1662 


1-V 


78 


Q 9+1-06 
^■■^_0 83 


0.29iSj'i^ 


1.47x10^° 


RXJ2302.8+0844 


0.722 


19 


918 


1-F 


108 




o-^7ttil 


4.91x10^° 


RXJ1113. 1-2615 


0.725 


19 


915 


1-F 


103 


6.411^2 
7.33+°1? 


0.35i°'f, 


5.48x10^° 


MS 1137.5+6624 


0.782 


20 


536 


1-V 


117 


0.24l°'i? 


1.2x10^° 


RXJ1350.0+6007 


0.804 


21 


2229 


1-V 


58 


4.381°-? 


0.57t°'^26 


1.78x10^° 


RXJ1317.4+2911 


0.805 


21 


2228 


1-V 


110.5 


3.771J-?? 




1.1x10^° 


RXJ17 16.4+6708 


0.813 


22 


548 


1-F 


51 




0-57lVit 


3.7x10^° 


MS 1054.4-0321 


0.832 


23 


512 


S-F 


80 




0.23lH^ 


3.6x10^° 


1WGAJ1226.9+3332 


0.89 


24 


3180 


I-V 


31.5 


12.03+^t? 


<0.11 


1.38x10^° 


CLJ1415. 1+3612 


1.013 


19 


4163 


I-V 


89 


0- 'O-0.67 

c or + 1.59 

J-OJ_1,08 


0.32i°;i« 


1.09x10^° 


RXJ09 10+5422 


1.106 


25 


2227, 2452 


1-V 


170 




2.2x10^° 


RXJ1252-2927 


1.235 


27 


4198, 4403 


1-V 


188.4 


6-81_o'87 


0.841^2? 


5.95x10^° 


RXJ0848.9+4452 


1.261 


26 


927, 1708 


1-V 


184.5 


A 44+1-11 

*-^^-0.86 


0.37+_lil 


2.6x10^° 


XMMUJ2235. 3-2557 


1.393 


28 


6975/6, 7367/8, 7404 


S-V 


196 


9.96_]^ 28 


OAltlit 


1.47x10^° 



Table 1 : The cluster sample used in this work, in order of increasing redshift. We show the properties of each data-set, includ- 
ing the references to previously published works, observation ID, detector (I=ACIS-I, S=ACIS-S) , observation mode (F=FAINT, 
V=VERYFAINT), and exposure time. The best-fit parameters T{zo) and Zpeizo) are obtained with the redshift frozen to the optical 
value Zo- Temperature and metal abundance refer to the global fit to the spectrum extracted within the radius Rext, and therefore 
represents an average value whenever a temperature or a metallicity gradient is present. The last column refers to the Galactic 
Hydrogen column density as measured by Kalberla et al. (2005). 

References. (l):Vikhlinin et al. (2005); (2):Martini et al. (2007); (3):Tucker et al. (1998); (4):Lewis et al. (1999); (5):Mushotzky & Loewenstein 
(1997); (6):Struble & Rood (1999); (7):Balestra et al. (2007); (8):Stocke et al. (1991); (9):B6hringer et al. (2000); (10):Ebeling et al. (2010); 
(ll):Moranetal. (2007); (12):Vikhlinin et al. (1998); (13):Kotov & Vikhlinin (2006); (14):Caccianiga et al. (2000); (15):Borgani & Guzzo 
(2001); (16):Mullisetal. (2003); (17):Allen et al. (2004); (18):Ebeling et al. (2007); (19):Perlman et al. (2002); (20):Gioia & Luppino (1994); 
(21):Holden et al. (2002); (22):Henry et al. (1997); (23):Tran et al. (2007); (24):Ebeling et al. (2001); (25):Stanford et al. (2002); (26):Rosati et al. 
(1999); (27):Rosati et al. (2004); (28):Mullis et al. (2005). 
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ID 


Zo 


zx 


T{zx) (keV) 




Cts (0.5-7 keV) 


SN 


Rext{") 




Abell907 


0.160 





1 cq+U. UU2 

— 002 







4g+U.04 

^"^ — 04 


48230±238 


202.5 


162 


4 


220.56 


Abelll689 


0.187 





1 QQ+0.005 
J^yy_0 005 


13.90+««^ 





Q 7+0. '07 
I _Q 07 


34008±200 


170.3 


93 


5 


31.59 


1E0657-56 


0.296 





qor+o.'ooe 


10 q+0.63 





9f^+0.'05 
^^ — Q 05 


35384±204 


173.1 


211 


6 


24.65 


MS1008. 1-1244 


0.306 





Q+0.014 
— 016 


6.05lgl 





97+0.'06 

06 


8993±103 


87.2 


108 


2 


24.32 


MS2137.3-2353 


0.313 





o-i n + 0.002 

0-i-0_Q QQ2 


5.26 12 





Qr+0:04 

oO_Q 03 


29509±176 


167.9 


76 


3 


171.03 


Abelll995 


0.319 





QI Q+0.003 

<JJ-0_Q QQ5 

QI r+0:005 
— 003 


9.051^3? 





41+0.07 

— 06 


27217±178 


152.9 


103 


3 


47.46 


ZW1358. 1+6245 


0.328 










^Q+0:07 
— 07 


18322±149 


123.2 


88 


6 


51.1 


MACSJ0404.6+1109 


0.355 


2 


^^^ — 048 




3 


09+if 


2368±59 


39.9 


120 


5 


3.67 


RXJ0027.6+2616 


0.367 





Q7Q+0.'024 
"J ' <J_0 023 







514^1! 


702±31 


22.7 


81 


2 


5.15 


MACSJ 1720.2+3536 


0.387 





QOC+O.OOS 
— 008 


6.5llg^ 





/lc+0.08 
^^_0 07 


6370±84 


75.9 


91 





50.46 


ZW0024.0+1652 


0.395 





QQ7+0.'006 

— 005 


4.58^Q27 





7(5+0.19 

— 17 


2671±59 


45.3 


59 





29.0 


RXJ1416.4+4446 


0.400 





^-1 Q+o:oo7 


c;+0.2'2 
—0 19 





o7+0.'25 
°'-0 21 


2120±51 


41.9 


76 


3 


31.91 


MACSJ0159.8-0849 


0.405 





4Q(5+o:oi 







qt;+0.08 
— 08 


7852±93 


84.5 


113 


2 


20.98 


MACSJ2228.5+2036 


0.412 





^-11 +0.008 

— 009 







qc+O.OO 
— OS 


6004±86 


69.4 


137 


8 


18.47 


MS0302.7+1658 


0.424 





^r,7+o:oi4 

1 —Q 014 


4.381^4 





45+0-24 
^'J — O 19 


652±27 


24.0 


68 


9 


8.48 


MACSJ0417.5-1154 


0.44 







12.581^^4 





QQ + O.'ll 


7544±95 


79.4 


132 


8 


9.84 


MACSJ 1206.2-0847 


0.44 












99+0.'08 
■^^_0 08 


11267±112 


100.2 


127 


9 


7.8 


RXJ1701. 3+6414 


0.453 





.c;4+0:008 
'-^ — 009 







r 1+0:13 
'-*J-_o 12 


2731 ±59 


46.1 


68 


9 


27.75 


RXJ 164 1.8+4001 


0.464 





+0.008 

'^'-*_0 01 







C/l+0:22 
— 19 


983±35 


28.4 


46 


7 


12.91 


MACSJ 1824.3+4309 


0.483 





368«-°» 


—0 96 





oq+o:48 

— 35 


469±26 


17.9 


78 


7 


7.9 


MACSJ131 1.0-03 11 


0.494 







8.79+S-o* 





QQ+0.15 

<jy_o 14 


2083±48 


43.6 


78 


7 


7.8 


RXJ1524.6+0957 


0.516 





[-90+0.013 

014 







04+0:14 

— 13 


1974±52 


37.7 


73 


8 


8.63 


4.6+0957 MS0015.9+1609 


0.541 





cr 1+0.009 

~0 008 


8.37+^'f 





Q9+0.'05 
— 05 


16826±139 


121.0 


118 


1 


49.81 


MACS J 1423. 8+2404 


0.543 





+0.012 

-0.013 


+0:71 

' • ' '-'—0.53 





01+0:1 

0-L_Q 09 


3518±62 


57.2 


76 


3 


14.07 


MACS J 1149.5+2223 


0.544 







12.75^Q gg 





27to-l 


9095±109 


83.3 


150 


1 


7.13 


SCI 120-1202 


0.562 





438l°-°l* 







C7+0.32 
^ ' —0 25 


714±33 


21.7 


51 


7 


6.48 


MS2053.7-0449 


0.583 





t-qc+o.oio 







9+0.11 
■^_0 11 


1884±49 


38.5 


49 


2 


3.69 


RXJ0956.0+4107 


0.587 





+0.013 
— 012 


5.79ll'g| 





7Q+b.81 

— 45 


467±24 


19.0 


64 





3.78 


MACSJ2129.4-0741 


0.589 





7+0.013 
^ ' —0 015 


WAtlil 





c7+0:i6 
^ ' —0 14 


2918±59 


49.1 


100 


9 


18.33 


MACSJ0647.7+7015 


0.591 





f^rr,+0;013 

^'-'■^ — O 013 







9.+o:i6 

19 


2899±58 


49.9 


78 


7 


3.52 


RXJ0542.8-4100 


0.634 





8.09lJ'^5 





iq+0.14 

— 13 


1908±50 


38.5 


64 





2.04 


MACSJ0744.9+3927 


0.698 





790+0.01 

— 013 


9.6llg-II 





Q^+o:o9 

— 08 


5669±79 


71.4 


81 


2 


17.68 


RXJ 122 1.4+49 18 


0.70 







S.ielfre 





Qc+o:i3 


2828±61 


46.6 


73 


8 


10.39 


RXJ2302. 8+0844 


0.722 





64tH'26 


6.6in'2 





9(5+0:i7 
^—0 15 


1387±44 


31.8 


46 


7 


3.27 


RXJl 113. 1-2615 


0.725 





77Q+b.014 
' ' — 015 


c 1+0.83 
—0 67 





t;c+o:21 
^^_0 18 


1135±38 


29.9 


39 


4 


13.42 


MS 1137.5+6624 


0.782 







7.46+gg 





04+0:12 


3957±67 


59.1 


44 


3 


10.59 


RXJ1350.0+6007 


0.804 





8341°-°?^ 


4.36lg-l3 





cc:+0.'36 
'-*'-*_0 27 


622±31 


20.3 


51 


7 


9.09 


RXJ1317.4+2911 


0.805 


1 


nKK+0.035 
— 031 


3.84lj'^^ 







230±20 


11.4 


29 


5 


3.05 


RXJ17 16.4+6708 


0.813 





00^+0.013 
009 









1341±40 


33.8 


44 


3 


19.28 


MS 1054.4-0321 


0.832 





Oq/,+0.021 


7 97+0.48 
—0 35 





90+0.07 

^■^-0.07 


8637±103 


83.8 


68 


9 


12.46 


WGA1226.9+3332 


0.89 





1Q7+0.027 
-0.032 


7 9r+0:85 
' •^'-'-0.69 





9+0.15 
^-0.13 


2353±51 


45.8 


66 


4 


2.55 


CLJ1415. 1+3612 


1.013 





Of. (-+0.034 
^"^-0.024 


6.311^:? 





00+0. 16 
■J"-0.14 


1275±39 


33.0 


36 


9 


8.53 


RXJ09 10+5422 


1.106 





/1SQ+0026 
^°''-0.022 


4 1+0.92 
^•^-0.56 





40+0.36 
^■^-0.25 


419±24 


17.2 


24 


6 


3.69 


RXJ1252-2927 


1.235 


1 


99f;+0.053 
^^"-0.03 


"•^-1.08 





04+0.32 
^4_o.28 


780±33 


23.4 


32 





14.56 


RXJ0849.9+4452 


1.261 


1 


or|4+0.054 


^••^-0.71 





r 1+0.48 
^J^-0.34 


327±22 


14.6 


22 


1 


2.87 


XMM2235. 3-2557 


1.393 


1 


000+0.046 
•J°'J-0.041 


Q (5+2.0 
^■"-1.12 





41+0.24 
^-■--0.21 


1433±45 


31.7 


24 


6 


4.25 



Table 2: Best-fit results for our sample of Chandra clusters obtained with a free redshift parameter The focus is on the best-fit X-ray 
redshift zx, compared to the optical value Zq- We also list the X-ray temperature T{zx) and the abundance Zpeizx)- Net photons 
(0.5-7 keV band) and S/N within the extraction radius Rext are also shown. The value of ACstat corresponding to the best-fit zx 
value is shown in the last column. 
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